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Abstract. We show that for the straightforward quantized rela- 
tivistic Coulomb Hamiltonian of a two-dimensional atom - or the 
corresponding magnetic quantum dot - the maximal number of 
electrons does not exceed twice the nuclear charge. The result is 
then generalized to the presence of external magnetic fields and 
atomic Hamiltonians. This is based on the positivity of 

|x|T(p)+T(p)|x| 

which - in two dimensions - is false for the non-relativistic case 
^(p) — P 2 /2, but is proven in this paper for T(p) = |p|, i.e., the 
ultra-relativistic kinetic energy. 

1. Introduction 

The energy of two-dimensional quantum systems interacting via three 
dimensional Coulomb-like potentials like graphene in a background 
magnetic field given by a vector potential 21 is given by the quadratic 
form 

(1) E[i>] = (xl>,W n i>) 

where ip is a sufficiently smooth element from the Hilbert space S)% 
built from states above the Fermi level: 

(2) ^^:=A + [L 2 (I 2 :C 2 )], A+ := X(o,oo)(A»). 

(Note that this model is physically relevant for large gaps and relatively 
small interactions. For strong interactions and small gaps a field theo- 
retic description allowing for particle-hole creation and annihilation is 
required (Paananen and Egger [10]).) The operator 

(3) W N := ^2 {T %n + ecr n ■ A(x n ) - e^(x n )) + ^ r 

n=l l<m<n.<7V ' X ' m X ™' 

is the multi-particle Weyl operator (massless Dirac operator) with ki- 
netic energy := a ■ (p + e2l) where — e is the charge of the electron. 
The Hamiltonian is defined via the quadratic form. The Hamiltonian 
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has been used to describe multi-particle effects of graphene quantum 
dots (see Egger et al pQ where its basic mathematical properties have 
been discussed). The electro- magnetic potentials (p and 21 are defining 
the quantum dot. 

To be concrete we mention that the choice used by Egger et al pQ (see 
also Paananen et al Qjp ) would be allowed, namely to take 21 = 2lo + A 
and A = with 



i.e., we have a homogeneous magnetic field of strength B orthogonal 
to the 2i-22-plane, and 



the potential that eliminates the magnetic field in a circle of radius R 
but leaves the field unchanged outside. In total this is a homogeneous 
magnetic field with a cavity of radius R around the origin. As electric 
field we could choose <p{x) = eZ/|x|, i.e., the potential of a defect atom 
placed at the origin. 

An alternative approach would be to take 21 = 2lo, i.e., define the 
vacuum with respect to the homogeneous magnetic field, and to pick 
A = A . Although for weak fields these two will yield similar results 
they are not identical. Nevertheless, by the variational principle for 
eigenvalues in gaps [3j, the latter will bound the former from below. 
In addition it allows for a more direct treatment of the excess charge 
problem. Because of this, we will direct our attention on the second 
choice. 

The corresponding quadratic form on /\ n=1 A + [C^°(IR 2 : C 2 )] is bounded 
from below, if and only if 



independently of the strength of the field 21 [?]. Egger et al [TJ studied 
also numerically how many electrons a quantum dot is able to bind in 
the context of a mean-field model. However the question of bound- 
ing the total number of electrons localized in the quantum dot was 
left unanswered. In the sequel we will address this question, how- 
ever not for the no-pair Hamiltonian defined by the quadratic form of 
Wn in A^Li #a but for the corresponding Chandrasekhar type opera- 
tor CA,tp,N, i-e., Wn with Ta replaced by |T&| self-adjointly realized in 



L 2 (R 3N : C 2 ) with domain H 1 ^™ : C 2 ). Again by the variational 



principle, the eigenvalues of this operator bound the eigenvalues of the 
no-pair operator from below, since the quadratic form of the no-pair 
operator is just a restriction of the Chandrasekhar case. 
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An a priori bound on the maximal number of electrons that can be 
bound has been derived by Lieb [S] for relativistic and non-relativistic 
Coulomb Hamiltonians in three dimensions. The idea of the proof is 
to multiply the Schrodinger equation by |x 1 |?/'(a;) and to integrate. It 
yields - in the atomic case - the bound N < 2Z + 1 and rests on the 
inequality 

(6) 0<|x|T + T|x|. 

For the non-relativistic kinetic T = p 2 /2 (jSJ) is equivalent to a Hardy 
inequality which is true in three dimensions but false in two dimension. 
The case of the relativistic kinetic energy T = |p| is reduced by Lieb 
[5] to the non-relativistic case. In other words, it is not clear a priori 
whether © holds. The purpose of this paper is to show this inequality 
and apply it to the Hamiltonian Ca,<p,n of the quantum dot. 
Before we actually do this we remark: 

• In three dimensions Nam [8] improved Lieb's result for the non- 
relativistic Schrodinger operator and Z > 6. He could show the 

i 

bounds iV < 1.22Z + 3Zz using similar ideas. 

• In two dimension, the above failure of positivity of flH]) can be 
controlled using an idea of Seiringer [T2] and Nam and Solovej 
[S] : Since the eigenvalues of the two-dimensional hydrogen atoms 
(two-dimensional Kepler problem) are — Z 2 /2(n + l/2) 2 (Fliigge 
and Marschall [21 Problem 24] each of multiplicity 2n+l (e = 1). 
This gives (ip, Ix^^ip) < 41og(Z 1 / 2 ) + 10 for the ground state 
of the N particle system which in turn yields 

(7) N <2Z + \og(Z 1/2 ) + ^. 

2. Main Inequality 

Our basic result is 

Theorem 1. Assume 21 G L 2 oc (IR 2 : C 2 ), p := -iV, andT % := |p + 2l| ; 
then 

(8) |x|T a + T a |x| > 
on C£°(R d ). 

Proof. By the diamagnetic inequality (Theorem [T81) 

(9) (^|p||0|)<^sgn(0)|p + 2l|0) 

we see that |x||p| + |p||x| > implies |x||p + 21 1 + |p + 2l||x| > 0. 
In other words, it suffices to prove positivity without the background 
magnetic field. Moreover: it suffices to consider > 0. 
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Next we recall the following representation of |p| (Lieb and Yau [5]) 
in position space 



(10) |p|<) = <* j dx / dy W<x)-V<y))(^)-,.( y) ) 

jRd jRd |X - yr +i 

with a d = r(^± 1 )/(27r^). By polarization of (HDD Inequality © is 
then equivalent to 



A- A* |x-y| J+1 

where we restrict to non-negative ■?/> because of the above remark. Now, 
setting ^ = g/\. \ d l 2 and regularizing to avoid the singularity at x = y 
we get 

, „ l*|ff(x) 2 i |y|g(y) 2 n (v\n(yj\ |y| + |x| 

t=lim»/ dx / dy- Wd |y|d ^ Jj(y, W^ 



|x - y| d+1 + 2^e(|x| d+1 + |y| d+1 ) 

:i:u / f 2 -'f'';l- |y| ' J/2 



+ - / dx / dy 



e-+o j Rd |x| d J R d | e -y|rf+i + 2— e(l + |y| d+1 ) 
1 

2_ 

> /■ dx ^)! r^(2- r i-^_ r -^ )(2r )-^ 

i« d l x l Jo r 

do; 



(1* 


i + |y 1) (</(*)- 


g(y)) 2 




x|<V 2 |x-y| 


d+ll 


y\d/2 



(12) 



>2-^ / dx^ r^(2r^-r^-r-^) 

./rod |x| ./q ^ 



do; 



1 ' — we 2 -f e (i 2 _|_ r 2 



2 



d+l 



c/(x) 2 Z" 00 dr d^i 



>2 2 / dx / — [r 2 +r 2 _ ( r 2 + r 2 )] 



r s >r- 

>0 



-1 \ d+l 1 1 %i d+i 



> 



— lot) 2 + e(- 2 +r2 



2 

t>3 



where e is any unit vector in Mr. The positivity of the bracket follows 
from the fact that the function f{a) := r a + r _Q is strictly monotone 
increasing for positive r. □ 

We remark that the above proof also shows that d = 2 is borderline 
for positivity. In fact, 

|x||p| + |p||x| > 2a d >y d 
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where 
(13) 

d-i f dr r ^1 _d = i i _i f dcu 
7rf = 2 ^ / _ r 2 +r 2 _( r2+r 2 / ^ 

which changes sign at d = 2 whereas |x||p| + |p||x| > 1 for d = 3. 



3. Application to 2d Quantum Dots 

In this section we consider a 2d quantum dot given by the Hamil- 
tonian 
(14) 

N e 2 
C a , v ,n ■= ^ [|p + 2l| n + ecr n ■ A(x„) - e^(x n )] + 2j 

n=l Km<n<N 



a simplified model of (j3J), self-adjointly realized in f)jv : = A n =i -^ 2 (^ 2 : 
C 2 ). Here 21 is the background magnetic field (j3j) and A is the mag- 
netic field defining the dot. Furthermore, we allow for an attractive 
essentially spherically symmetric attractive potential. We have 

Theorem 2. Assume 21 E Lf oc (R 2 : IR 2 ) and |A(x)| < e<J/|x|, <p(x) < 
eZ/\x\, Z G [0, Kf.]. Assume that Ca,<p,n has a ground state with ground 
state energy E N below the saturation threshold, i.e., E N < E N _ X . Then 

N <2(5 + Z) + l. 

Note that our bound - in the absence of an electric potential - grows 
linearly in the missing magnetic flux fi which, in the case of (jSJ), equals 
fi = e 5 = R 2 Be/2. 

Equipped with Inequality (jSj) the proof follows now the lines of Lieb's 
original proof. 

Proof. Assume that there is ground state ip of Ca, v ,at, i-e., 
(15) C A ,<p,Nip = E N ip 

with En- We begin by singling out the first coordinate and multiplying 
(TT5|) by |xi|^> and obtain using dHJ) and a standard density argument 
(approximating if) by functions) 



(16) e 2 (|x 1 |^,[-(5 + Z)/|x 1 | + ^^^]^) 

< [Ca, v ,i - e 2 l x i _ XnT 1 ]^ + 1 ® (Ca^n-i - E N „ x )if)) 

= (^-^-Odxil^.V) < o. 
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Repeating the same argument for all the other coordinates x 2 ,...,x r 
and summing the result, gives 

1 N 

(17) -(5 + Z)N + -J2 



2 

m,n — 1 



Thus the triangle inequality implies 

-(5 + Z)N < 

which is the claimed result. □ 
Appendix A. Auxiliary Results 





+ x n | 


l^m 




N 2 - 


N 


2 



< 0. 



Theorem 3. For<j),r] e H 1 (E <i ), rj > 0, m e R+, T m (p) = a/p 2 + m 2 - 
m, 21 G L 2 3C (M 3 ) and p = — iV we have 

(18) {r,,T m (p)\<f>\) < ^,sgn(0)T m (p + 2t)0). 

This is Kato's inequality for the Chandrasekhar operator (Simon 
[T3] . Konenberg et al [I]). 
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